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The Flory \ parameter is typically assumed to depend only on the temperature, T . Experimental 
results often require the replacement of this x(T) by Xe-ffi that depends also on the monomer 
volume fraction, (j>, Xef f(<l>,T). Such Xef f(4>^T) can arise from two state models, proposed for 
polyetheleneoxide (PEO) and other neutral water-soluble polymers. The predicted (j> dependence of 
X = Xeff ~ (1 ~~ 0)^Xe/ / /d<f>, obtainable from colligative properties, differs qualitatively between 
the various models: (i) The model of Karlstrom (J. Phys. Ghera. 1985, 89, 4962) yields dx/d(f> > 
while the model of Matsuyama and Tanaka (Phys. Rev. Lett. 1990, 65, 341) and of Bekiranov et 
al (Phys. Rev. E 1997, 55, 577) allows for dx/d(j> < (ii) x((f>) as calculated from the Karlstrom 
model, utilizing the parameters used to fit the phase diagram of PEO, agrees semiquantitatively with 
the experimental values. On the other hand, x(0) similarly calculated from the model of Bekiranov 
et al differs qualitatively from the measured results. Altogether, %(</>) provides useful measure for 
the performance of a model. 



I. INTRODUCTION 

The Flory free energy plays a central role in the ther- 
modynamics of polymers. In this free energy, the mixing 
energy term has the form xfiO- — 4>) where \ is the Flory 
interaction parameter and </> is the monomer volume frac- 
tion. Typically, it is assumed that \ depends only on 
the temperature, \ — v(T).(l], ||, |], Q However when 
thermodynamic data are analyzed in terms of the Flory 
free energy, it is often necessary to replace x&O- ~ 4 1 ) by 
Xe//<^(1 — <t>) where Xeff is a function of both T and 
<f> i.e., xeff = Xeff{T,(j>)% I §, 0, § Considerable 
effort was devoted to clarifying the statistical mechan- 
ical origins of this and other deviations from the Flory- 
Huggins theory. Various extensions of the Flory-Huggins 
theory allow for compressibility, p local composition! 10 
and the role of intrachain contacts. [jll| Recently, the lat 
tice cluster theory utilized a more accurate solution of 
the lattice model while allowing for the structure of the 
monomers. |l2| The Flory-Huggins theory was also ex- 
tended in order to account for the phase behavior of 
aqueous solutions of neutral water-soluble polymers and 
in particular polyetheleneoxide (PEO). In such systems 
the solubility of the polymers is attributed to the forma- 
tion of H-bonds. These extensions introduce a two-state 
model wherein the monomers of the chain can exist in 
two distinct and interconverting states (l^, [l4|, [0| [Ig| 
(Figure 1). In the following we demonstrate that "two- 
state" models yield Xeffifi) and analyze the </> depen- 
dence of the resulting Xeff- Our discussion is mostly 
concerned with two models: The one proposed by Karl- 
strom (K model) (l^l and the model of Matsuyama and 
Tanaka Jl5| as formulated by Bekiranov, Bruinsma and 
Pincus (MB model). (T| 

The two models differ, as we shall discuss, in the pre- 
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Figure 1. A schematic representation of the dynamic 
equilibrium in two-state polymers where the monomers exist 
in two interconverting states (inset), A (empty circle) and B 
(full circle). 

cise identification of the monomeric states. However, in 
both cases it is assumed that one monomeric state, A, is 
"hydrophilic" while the other, B, is "hydrophobic". As a 
result, a binary solution of a polymer in water behaves in 
fact as a reactive ternary system. In particular, the frac- 
tion of A monomers, P, depends on both T and 4>. Con- 
sequently the different versions of this model yield, as we 
shall see, Xef / that depends on both T and 4>. The orig- 
inal papers proposing theses models |l4|, [l^, focused 
on the phase behavior while overlooking the emergence 

of X e//(0)@. 

The goal of our work is to obtain explicit expressions 
for Xeff{4>) f° r the K and MB models, analyze the <f> de- 
pendence and compare the results to experimental data. 
It is noteworthy that by allowing for the internal degrees 
of freedom of the monomers, the two-state models can ra- 
tionalize the observed Xe//(</>) while retaining the main 
approximations of the Flory-Huggins theory i.e., incom- 
pressibility, monomers and solvent molecules of identical 
size and shape and local concentration that equals the 
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Figure 2. Plots of measured \ — Xeff ~ (1 — 0)^Xe/ /d<j> vs. 

<j> for (a) aqueous solutions of the neutral water soluble 
polymers PEO N = 3000, T = 338K (full circle, afterH |) 

and PVP T = 295K (empty circle, after[[t|]) and for (b) 
nonaqueous solutions of polyisobutylene in benzene (dashes) 
and of polystyrene in toluene (continuous line) pi . 

global one. Aside from the fundamental interest in the 
origin of the <f> dependence of Xeff this analysis under- 
lines the utility of Xe//(0) as a criterion for the perfor- 
mance of the models. 

While both the K model and the MB model recover the 
observed phase diagram of PEO, only the K model yields 
a qualitatively correct Xeff '{<f>)- O n the other hand, an 
opposite trend is found for Xeff{4>) of polyvinylpiridon 
(PVP) which exhibits a <\> dependence that differs qual- 
itatively from that of PEO |]l9) (Figure 2). In this case 
the K model is inherently incapable of reproducing the 
experimental results while the MB model does allow for 
a similar behavior. 

Our analysis is mostly concerned with the concentra- 
tion dependence of 



X = Xeff - (1 



dx 



eff 



(1) 



We choose to focus on x because measurements of col- 
ligative properties, such as vapor pressure and osmotic 
pressure, yield x rather than Xe//-B As a result, it is 
the x values that are usually reported in the literature. 
Furthermore, it turns out that the <p dependence of x(0) 



is simpler to analyze. Since the different versions of the 
two-state model aimed to rationalize the phase behav- 
ior of aqueous solutions of PEO, the parameters involved 
were chosen with view of recovering the known phase di- 
agram. However, these models are in fact suitable candi- 
dates for the description of aqueous solutions of neutral 
water-soluble polymers in general. With this in mind 
we consider the two-state model as a generic model for 
such systems rather than as a specific model for PEO. 
Accordingly, we make no specific assumptions regarding 
the choice of the various parameters involved. 

The paper is organized as follows. A unified descrip- 
tion of the K and MB models, stressing the similarities 
and the differences, is presented in section II. Section 
III contains a brief summary of the thermodynamics of 
polymer solutions with Xeffiffi)- This discussion intro- 
duces the approach utilized later for the study of xi'P) ■ 
Xeff(<ft), x(4") an( i dx{4>)/d(t> for the K and MB models 
are discussed in sections IV and V respectively. In the 
discussion we compare the K and MB models with two 
other physically transparent models yielding a <p depen- 
dent Xeff- 

II. TWO "TWO-STATE" MODELS 

In the "two-state" models for water-soluble homopoly- 
mers the monomers of the chain may assume two inter- 
converting forms characterized by different interaction 
energies. In particular, one monomeric state, A, is as- 
sumed to be "hydrophilic" while the other, B. is assumed 
to be "hydrophobic" . The K and MB models differ some- 
what with respect to the underlying physics. The K 
model HH focuses on molecular rotations giving rise to 
monomeric states of different polarities. The state en- 
dowed with a strong dipole moment is identified as the 
hydrophilic species, A, while the less polar state is identi- 
fied as the hydrophobic state, B. Within this model the 
interconversion is a unimolecular reaction, A £?, and 
the mass action law is 
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Here P is the fraction of monomers in a hydrophilic A 
state and Kk is the corresponding equilibrium constant. 
The explicit form of Kk, to be given in section IV, is irrel- 
evant at this point. The MB model|l7j focuses on the role 
of H-bonds. In particular, a monomer that forms H-bond 
with water is considered, in effect, as hydrophilic while a 
monomer that does not form such bond is hydrophobic. 
Within this model one distinguishes between free water 
molecules (S) and water molecules that are bound to the 
polymer chain. This distinction has no counterpart in 
the K model. Accordingly the interconversion reaction 
is a bimolecular reaction, B + S A, and the chemical 
equilibrium is specified by 
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where <po is the volume fraction of water. The explicit 
expression for Kmb, to be given in section V, is not re- 
quired for the present discussion. Thus, in the K model 
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while in the MB model 



= 1 



P(j) 



(4) 



(5) 



where P is the fraction of monomers that form H-bonds 
and Pip is the volume fraction of bound water. As a 
result, the two models involve different reference states. 
In particular, the melt state, <j) = 1, in the K model is 
characterized by 



P{cj> = 1) = P* > 
while in the MB model 
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Both models utilize a Flory type approach to obtain the 
thermodynamics of the solution. In particular, random 
mixing is assumed and the energy per site of a polymer- 
water solution inscribed on a lattice is written as 
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where z is the coordination number of the lattice and 
k is the Boltzmann constant, caa> ebb an d ess denote 
respectively the energies of A, B and solvent molecules 
in a pure bulk phase, xab, Xas and xbs denote the 
Flory interaction parameters corresponding respectively 
to AB, AS and BS interactions (where xas — £as — 
{f-AA + ess)/2 etc). The first term allows for the energy 
of the A and B in a hypothetical pure bulk phase while 
the second reflects the energy of the solvent molecules 
in a pure bulk state. The interactions between the two 
monomeric states and the solvent gives rise to the third 
term while the last term allows for the AB interactions. 
For future reference it is convenient to define a dimen- 
sionless parameter 



Ae = 



£AA - £BB 

' 2kf 



(9) 



characterizing the difference between the energies of the 
two monomeric states. This expression reduces to the 
familiar Flory form when P = 1 or P — that is, when 
the monomers have only one state. The energy E is sup- 
plemented by two entropy contributions per lattice site. 
One term allows for the different possible sequences of A 
and B states along the chain 



Sab 



>[PlnP+(l-P)ln(l-P)] (10) 



while the second reflects the translational entropy of the 
chains 
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The Sab term is distinctive to the two-state models. It 
vanishes when the monomers exist in a single state i.e., 
P = 1 or P = 0. As expected, Straus is identical in form 
to the entropy term in the Flory free energy. Altogether, 
the free energy of the solution is F = E — T (Sab + Straus) 
and the mixing free energy is 

F mix = F(<j>, P) - 4>F{4> = 1) - (1 - 4>)F(4> = 0) (12) 

Both models neglect the fluctuations in the number of 
A monomers per chain, m. Strictly speaking, chains with 
different m are distinguishable and should be treated 
as different chemical species. [[L5[ Rather than allow for 
this multiplies equilibrium the K and MB models con- 
sider a single polymer species characterized by an av- 
erage m — PN where iV is the polymerization degree. 
This approximation is justified because the m values are 
described by a sharply peaked Gaussian distribution. Jljj] 
While the two models lead to similar free energies, the 
free energies invoked in the K and in the MB-models 
differ in two respects. In the MB model (i) The trans- 
lational entropy is P dependent because of the <po In </>o 
term and, (ii) since P(<j> = 1) = 0, F m ix does not reflect 
contributions due to P(<f> = 1) = P* > that appears in 
the K model. 

Finally, it is helpful to note certain additional differ- 
ences between the models as described in the original 
papers. |lj, [TtJ While the physical content of these is mi- 
nor in comparison to the differences discussed above, they 
are of interest for comparison purposes. Thus, the anal- 
ysis of Bekiranov et al focused on the particular case of 
XAB = 0, xas = Xbs = X = 2.885 - 0.0036P and Ae = 
6.38-2408/T while in the K model xas = 80.0/T, xbs = 
684.5/P, xab = 155.6/T and Ae = -625.2/T + In 8. 
The temperature dependence of the x parameters in the 
two papers is somewhat different. Bekiranov et al as- 
sume x ~ a + bT while in the K model x ~ b/T. In 
both models Ae ~ a + b/T. In the model of Bekira- 
nov et al this follows from the identification of Ae with 
AF /kT where AP = AP - TAS is the free energy 
change associated with the formation of H-bond. The 
a + b/T dependence in the K model is a consequence of 
the degeneracy of the monomeric states. It assumed that 
each of the two monomeric forms can exist in a number 
of equivalent states and the resulting entropy gives rise 
to a. The analysis of Bekiranov et al allows for the ef- 
fect of pressure. This is attributed to the existence of 
a preferred volume for a H-bond. 63] The application of 
pressure is assumed to introduces a geometric constraint 
resulting in a lower number of bound water molecules. 
This argument yields a semiquantitative agreement with 
experimental results on the phase diagram of PEO. Note 
that this coupling mechanism can only operate within the 
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MB model and can not be used in the K model. Further- 
more, it is not fully consistent with the incompressibility 
assumption invoked by Bekiranov et al. 



III. THERMODYNAMICS OF POLYMER 
SOLUTIONS WITH Xeff(4>) 

Before we proceed to analyze the <f> dependence of x 
within the two "two-state" models, it is helpful to sum- 
marize the thermodynamic relationships that apply when 
X = x(T) is replaced by X eff = Xeff(T, The Flory 

like mixing free energy per site, F mix , consists then of 
two terms. One is an interaction free energy that is the 
counterpart of the mixing energy x4>{^ ~ 4>) 



F- 

1 1:1 



kT 



Xeff WW -<!>)■ 



The second is the translational free energy 
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The chemical potential of the solvent is u s = /x° (P, T) 
7ra 3 while the osmotic pressure 7ra 3 = 4> 2 A ( F "' 3! ) is 



7ra 
~kT 



N 



ln(l 



-x<f 



(15) 



de- 



where x = Xeff — (1 — 4 1 ) dX d<p f replaces x- Since fi s 
termines the colligative properties of the solution, mea- 
surements of such properties yield x rather than Xeff - I n 
turn, x is obtainable from Fi nt via 
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d F, 
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In the next two sections we will utilize equation ( |l6| ) 
to obtain x from F int = F mix - F tra ns- This equation 
also shows that terms of the form const' <p in Fi nt do not 
contribute to x- This will allow us to ignore such linear 
terms that arise because of the choice of the reference 
state. 



IV. x WITHIN THE K MODEL 



Within the K model the mixing energy per lattice site, 

Fmix j is 



E,, 



kT 



- 0PAe + 0(1 - 4>)[P X AS + (1 - P)xbs] + 
4> 2 XabP(1 — P) + terms linear in (j) (17) 



The first term measures the difference in the energy of 
isolated A and B monomers. The second term accounts 
for the interactions between the monomers and the sol- 
vent. It is a generalization of the x4>(l ~ 0) term in the 




Figure 3. \k vs - 4> (continuous line) calculated using 
X as = 80.0/T, X bs = 684.5/T, X ab = 155.6/T and 
Ae = — 625.2/T + In 8, the parameters used to fit the phase 
diagram of PEO, agree semiquantitatively with the 
experimental values of x (full circles) as measured with 
N = 3000 at (a) T = 328K and (b) T = 338-ff. While the 
numerical values and the curvature of Xk change with the 
choice of xas, Xbs, Xab and Ae, it is normally a 
monotonically increasing function of <f>. 



Flory free energy for a binary solutions. The third term 
allows for the interaction between A and B monomers. 
The mixing entropy is 
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= -4>[P InP + (1 - P) ln(l - P)} - 
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terms linear in 



The first term allows for the AB mixing within the chains 
and the second for the translational entropy of the chains. 
In both expressions, the terms linear in (ft arises because 
the free energy of the melt state depends on P* and 
F m ix = F - (1 - <t>)F{<t> = 0) - 0P(0 = 1,P»). Upon 
discarding terms linear in 0, F int = F m %x - F tra ns is 
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The equilibrium value of P for a given </> and T, P e = 
P e {<j>), is specified by the condition dF int /dP = leading 
to 



P 



= K 



K\ 



(20) 



1-Pe 

exp [-Ae - (1 - 4>){xas - Xbs) ~ ^XAb(1 - 2P e )] 

where Kk is the equilibrium constant for the unimolec- 
ular intrachain A B reaction. When \as — Xbs = X 
and xab = 0, the parameters chosen by Bekiranov et. 
al, Kk is independent of (j). Consequently, for this case 
P e is independent of </>, dP e /dcf) = 0, thus leading, as 
we shall see, to dx K /d(j> = 0. To assure the stability of 
the A B equilibrium we further invoke the condition 
8 2 F mt /dP 2 > 



d 2 F mi 
dP 2 
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P(l - P) 
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(21) 



This ensures that the free energy curve is concave and 
that fluctuation relax back to the equilibrium state ac- 
cording to the LeChatelier principle. Using x— ~~§$vr$ 
and the equilibrium condition we find 

XK = PeXAS + (1 - Pe)XBS ~ XA B P e (l ~ Pe) (22) 

Thus, within the K model the <fi dependence of Xk is 
due to P e (</>)■ In order to specify dx K /d<t> it is first nec- 
essary to obtain dP e ((f>) /d(j> by differentiating the equi- 
librium condition (|20|) with respect to <fi leading to 



dP e {<j>) _ SfXK 

dcj) 



dP e d 2 F mt /dP 2 



(23) 



Here, and later in (p^), d 2 Fi nt /dP 2 signifies the value 
of (|2lJ) in equilibrium i.e., P = P e {4>). Equation (§|) 
indicates that the sign of dP e (4>) / 'dcj) is determined by 
dx K /dP e . This reflects the concentration dependence of 
Kk = const' exp((f>dx K / dP e ) . Altogether 



dx K _ dXK dp e 



dP K 



K 



dP e J d 2 F int /dP 2 



> (24) 



and Xk is an increasing function 4> for any choice of the 
parameters xas, Xbs, Xab and Ae except for the case 
Xas = Xbs and xab = when dx K /d(j) = 0. The 
calculated Xk W j utilizing the parameters used to fit the 
phase diagram of PEO, is in semiquantitative agreement 
with the experimental results for %(</>) (Figure 3). The 
resulting curve is representative of the behavior of Xk {4>) 
in general. 



V. x WITHIN THE MB MODEL 

The mixing energy per site in the MB model is essen- 
tially identical to that of the K model 
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Figure 4. Xmb vs - plots exhibit qualitatively different 

forms, depending on the choice of parameters (a) 
monotonically increasing xas = 0.7, Xbs = 0.9, xab = 2 

and Ae = 4, (b) monotonically decreasing xas = 0.1, 
Xbs = 0.4, xab = -3 and Ae = 5, (c) maximum \as = 9, 
Xbs — 0.1, xab — 0.1 and Ae = —5, (d) minimum xas = 3, 
Xbs = 3, xab = 0.5 and Ae = -2. 
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= 0PAe + (1 - 0)0 [P X as + (1 - P)xbs] 



4> 2 XabP{1 - P) 
while the mixing entropy per site is 

: -0[PlnP+(l-P)ln(l-P)] 



(25) 
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In cf) + (1 - (j) - Pcj>) ln(l -<j)-P<j)) 



(26) 



These differ from the corresponding expressions within 
the K model in two respects: (i) Terms linear in (j) are 
no longer present because the dependence on P* disap- 
pears since P((j> — 1) = (ii) The translational en- 
tropy contribution to S m i x now depends on P because 



= 1 - 



P0. Fi, 



Ftrans IS 



kT 



= <j>PAe + <t>(l-<f>)[P X AS + (l-P)XBs} + 

2 XabP(1 -P)+ 0[PlnP + (1 - P) ln(l - P)] - 

P0 



P01n(l-0) + (l-P0-0)ln 1- 



1 - 



(27) 



The equilibrium value of P, for a given <j) and T is 
specified by the condition dFi nt / dP = 



K 



AIB{ 



>>,Pe,T) 



P 



(28) 



(l-P e )(l-P e 0-0) 

exp [1 - Ae - (1 - 4>)(xas ~ Xbs) - ^Xab{1 - 2P e )] 

where Kmb = cKk is the equilibrium constant for the 
B + S A reaction. As in the K model, Kmb is inde- 
pendent of when xas = Xbs and xab = 0. 

However, in the MB model P e depends on cj) even when 
Km b does not because of the loss of translational entropy 
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Figure 5. Xmb vs - 4* (continuous line) calculated using 
Xas = Xbs = X = 2-885 - 0.0036T, X ab = and 
Ae = 6.38 — 2408/T, the parameters used to fit the phase 
diagram of PEO, qualitatively differ from the experimental 
values of x (full circles) as measured with N = 3000 at (a) 
T = 328K and (b) T = 338K. 

of the bound water. The stability condition for the B + 
S ^ A reaction, d 2 F int /dP 2 > 0, leads to 



d 2 R 



iiit 



dP 2 

Using x 
we find 



P(l - P) 



2<t> X 



AB 



1 - 



P(j) 



> 0. 
(29) 

|#r and the equilibrium condition (|2q) 



d<f> kT<p 



Xmb = PeXAs + (I ~ P e )XBS -XabPc{1 - Pe) 

Pe , 1 . / P e 4 



(30) 



While the first three terms are identical in form to Xk 
these are now supplemented by two additional terms aris- 
ing from the translational entropy of the water. When 
XAS = Xbs = X and Xab = the interaction terms are 
constant, Xk — X, and the <fi dependence of Xmb re ~ 
fleets solely the entropic contribution. Again we obtain 
dP e (4>)/d(j) by differentiating the equilibrium condition, 



|), with respect to <j> finding 
dP e {4>) 9xmb 



dd 



dP e d 2 F mt /dP 2 



(31) 



In equation @, and later in (||), d 2 F int /dP 2 de- 



notes the value of ( |29| ) in equilibrium. While equa- 
tion (|l]) is similar in form to (^3|) the two are not 
identical because the two models yield different P e and 
because 9% MB /9P e ^ &x K /dP e . As before the sign 
of dP e (4>)/d(j> is determined by dx MB /dP e - This re- 
flects the concentration dependence of K MB = (1 — 
4> — P4>)Kmb in particular, dx MB /dP e = d\a.K MB /d4>. 
Upon substituting dP e (4>)/d<fi into dx MB /d4> as found 
from PQ) we obtain 



dx_ 



MB 



dx 



MB 



dP, 



d 2 F mt /dP 2 
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1 



(l-0)(l-P e 0-0) 



(32) 



The first two terms are positive. Of these, the first 
is analogous to the corresponding result within the K 
model. The third term is negative. As a result, 
dx MB /d(j> is no longer positive definite. Depending on 
the choice of xas, Xbs, Xab and Ae the calculated Xmb 
can be monotonically increasing or decreasing as well as 
exhibit a minimum or a maximum (Figure 4) . In the case 
of PEO the calculated Xmb (0) > utilizing the parameters 
used to fit the phase diagram of PEO, differs qualitatively 
from the experimentally measured x{4>) (Figure 5). 



VI. DISCUSSION 

As we have seen, two versions of the "hydrophobic- 
hydrophilic" two-state model for water-soluble polymers 
in aqueous solutions lead to a <fi dependent Xef / • The In- 
version, where the two states undergo unimolecular intra- 
chain conversion, results in dx/dcj) > 0. Using xas, Xbs, 
Xab and Ae as obtained by fitting the phase diagram of 
PEO this model yields x that agrees semiquantitatively 
with the experimentally observed values || Q (Figure 3). 
Within the MB model the hydrophilic monomeric state 
binds a water molecule. As a result, the interconver- 
sion reaction is bimolecular and the translational entropy 
of the water plays a role in determining the equilibrium 
state. In this model x can display a number of scenarios 
(Figure 4): x(0) ca n be monotonically increasing, mono- 
tonically decreasing or exhibit an extremum (maximum 
or minimum). 

The x vs. 4> as calculated with the parameters used to 
fit the phase diagram of PEO, differs qualitatively from 
the experimentally obtained curve(Figure 5). However, 
the dx/d(f> < behavior allowed by this model is of inter- 
est since it has been observed in aqueous solutions of the 
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neutral water-soluble polymer PVP[fL9| (Figure 2). These 
results stress the importance of using the experimental 
values of x in fitting xas, Xbs, Xab and Ae and in eval- 
uating the performance of the model. From the perspec- 
tive of the general theory of polymers it is of interest that 
the "two-state" models can account for the dependence of 
Xeff on 0, T and the pressure while retaining the prin- 
cipal approximations of the Flory-Huggins theory i.e., 
incompressibility, monomer and solvent of identical size 
and shape and local composition that equals the global 
one. 

At this point it is helpful to consider alternative mech- 
anisms leading to Xeff (0)- For brevity we limit this 
discussion to two physically transparent models. One 
is the n-cluster model proposed by de Gennes for aque- 
ous solutions of PEO.[jl6| This model was motivated by 
two experimental observations concerning the behavior 
of such solutions: (i) the interpretation of calorimetric 
measurements in terms of the Flory mixing free energy 
yields a % that increases strongly with 0,|| and (ii) re- 
ports of formation of aggregates in concentrated solu- 
tions of PEO.[^l| Within this model the concentration 
dependence of Xeff is attributed to attractive interac- 
tions leading to stable clusters of n > 2 monomers while 
binary monomer- monomer interactions remain repulsive. 
This is another variation of a two-state model involving 
a dynamic equilibrium between "clustered" and "unclus- 
tered" monomers. In molecular terms, a n-cluster may 
correspond to a micelle or a mixed helix. The forma- 
tion of the n-clusters gives rise to an additional term, 
—p(T)cj) n (with p > 0), in the interaction free energy. 
Thus, the mixing free energy per site is 



^ n 



kT 



x 0(l-£)+p(T)(0-0") + 



^ln^+ (1 - 0)ln(l 



leading to 



Xeff = X + P- 



1 - 



The corresponding x is 

X = X + p{n- IW 1 ' 2 



(33) 



(34) 



(35) 



and dx/dcf) > 0. The n-cluster model is indeed capable 
of rationalizing the two experimental observations noted 
above. However of these, the second observation is now a 
subject of debate. |22[ Furthermore the molecular struc- 
ture of the ethylenoxide monomers does not reveal am- 
phiphilic motifs. It is thus difficult to justify the as- 
sumption of cluster formation in solutions of PEO. Ac- 
cordingly, the validity of this model in the case of aqueous 
solutions of neutral water soluble polymers such as PEO 
is not obvious. On the other hand, the model is indeed 
applicable to solutions of polysoaps where formation of 
inter and intrachain micelles does occur |E3|. The second 



model, advanced by Painter et al, is applicable to all poly- 
mer solutions. Within this model, the dependence 
of Xeff is attributed to the interplay of intrachain and 
interchain contacts. The authors argue that there is a 
probability 7 for intrachain monomer- monomer contacts. 
As a rough approximation, 7 can be identified with the 
monomeric volume fraction, cpQ, within a Gaussian coil 
of radius Rg ~ N 1 / 2 a where N is the polymerization de- 
gree and a is the monomer size, 0g ~ Na 3 /R^ « N^ 1 / 2 . 
In a lattice comprising of Nt sites of coordination num- 
ber z there are Nj-cj) sites occupied by monomers with 
a total of Ntz4> adjacent sites. Of the Nxzcf) adja- 
cent sites Nxz(j)"f are occupied by monomers because 
of intrachain contacts. The total number of free, un- 
blocked adjacent sites is thus Ntz(1 — <jry) while the num- 
ber of free sites adjacent to monomers is Nxzcf)(l — 7). 
The probability of a free site adjacent to a monomer 
is thus 0(1 — 7)/(l — 07). Accordingly, the number of 
monomer-solvent contacts is Nt(1 — 0)0(1 —7)/(l — 07). 
This expression allows for the requirement that a sol- 
vent molecule occupies an unblocked site adjacent to a 
monomer. Accordingly, the mixing energy term per site 
is E mlx /kT = x(l ~ 0)0(1 ~ 7)/(! ~ 07) yielding an 
enthalpic contribution to Xef f 



1-7 

and to a corresponding enthalpic x 



X = X 



1-7 
1-07 



(36) 



(37) 



where, again, dx/dcf) > 0. In their original paper, Painter 
et al supplemented this enthalpic x by an entropic one, 
allowing for the effect of the chain bending back on itself 
following the analysis of Huggins.[Q 

Aside from physical insight regarding the molecular 
origins of Xeffi'f), the above discussion identifies cer- 
tain difficulties. Overall, it seems unlikely that one of 
the four models described above will emerge fully victo- 
rious. By construction, the n-cluster model is applica- 
ble only to solutions of associating polymers. The two 
"two-state" models are suitable candidates for the de- 
scription of aqueous solutions of neutral water-soluble 
polymers exhibiting insolubility gap. While the model of 
Painter et al applies, in principle, to all polymeric systems 
it can not account for systems exhibiting dx/dcf) < 0. 
Such behavior was actually observed in both aqueous 
and non-aqueous solutions, yet of the four models consid- 
ered above only the MB model yields a scenario involv- 
ing dx/dcf) < 0. Altogether one can thus envision situ- 
ations where all the different mechanisms described may 
contribute simultaneously. Another discouraging obser- 
vation concerns the number of parameters involved. All 
four models introduce additional parameters that do not 
appear in the familiar Flory-Huggins theory: 7 in the 
model of Painter et al, p and n in the n-cluster model, 
Ae, xab, Xas and xbs for the two "two-state" models. 
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The necessity to unambiguously determine the additional 
parameters limits the predictive power of the models. 
Clearly, this problem is even more serious when a number 
of mechanisms contribute simultaneously to the <j> depen- 
dence of Xef f ■ Note however that this last difficulty can 
be partially resolved in certain cases. Thus the contribu- 
tion of the mechanism of Painter et al can be separated 
from the one due to the "two-state" mechanism. This is 
because the mechanism of Painter et al is inherently a 
polymeric effect that disappears in the monomeric limit, 
N = 1. On the other hand, the two-state models also ap- 
ply to solution of the unpolymerized monomers i.e., Xeff 
and x are independent of N. 

Following our discussion of alternative models for 
Xeff ((ft) it is helpful to summarize the evidence favoring 
the two-state models. Using these models it was possi- 
ble to derive the phase diagram of PEOjlJ, |l7j and its 
pressure dependence. |l7|] In addition, the K model was 
successfully utilized to rationalize a variety of systems in- 
volving PEO, Polypropyleneoxyde (PPO)and PEO-PPO 
copolymers. 1 25 These successes lend support to the va- 
lidity of the two-state model. On the negative side, 
there is relatively little direct evidence for the exis- 



tence of the two monomeric states. Computational stud- 
ies lend the most direct support to the two-state mod- 
els. The K model is supported by recent molecular 
dynamics study that monitored the population of hy- 
drophilic and hydrophobic conformations as a function 
of temperature and concentration and relates their rela- 
tive stability to polar interactions. |2(| This study utilized 
force fields obtained from quantum mechanical calcula- 
tions of dimethoxyethane. [g7J Somewhat different con- 
clusions were reached in studies of higher oligo (ethylene 
oxides) . Q In particular, these studies (i) support the 
MB model (ii) stress the importance of conformers al- 
lowing formation of "water bridges" i.e., where a single 
water molecule forms two hydrogen bonds with the oxy- 
gens in the oligo (ethylene oxides). One such conformer, 
though not the most stable, is helical, (iii) allow to ratio- 
nalize force measurements of PEO[^9| and, in particular, 
the difference in the extension force laws found in water 
and in non- aqueous solvent. [pT)| Thus while a consen- 
sus on the computational results is yet to emerge, this 
approach may eventually allow to identify the relevant 
monomeric states and to obtain Ae, xab, Xas an( A Xbs- 
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